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Abstract 

We apply the semiclassical spin coherent state method for the den- 
sity of states by Pletyukhov et al. (2002) in the weak spin-orbit cou- 
pling limit and recover the modulation factor in the semiclassical trace 
formula found by Bolte and Keppeler (1998, 1999). 



1 Introduction 

A new solution to the problem of how to include spin-orbit interaction 
in the semiclassical theory was recently proposed by Pletyukhov et al. 
Q . They use the spin coherent states to describe the spin degrees of 
freedom of the system. Then a path integral that combines the spin 
and orbital variables can be constructed, leading to the semiclassical 
propagator (or its trace) when evaluated within the stationary phase 
approximation. In such an approach the spin and orbital degrees of 
freedom are treated on equal footings. In particular, one can think 
of a classical trajectory of the system in the extended phase space, 
i.e., the phase space with two extra dimensions due to spin. (The 
spin part of the extended phase space can be mapped onto the Bloch 
sphere.) Like in the pure orbital systems, it is possible to construct 
a classical Hamiltonian that will be a function of the phase space 
coordinates. The trajectories of the system satisfy the equations of 
motion generated by this Hamiltonian. 

In this letter we apply the general theory Jjj] to the limiting case of 
weak spin-orbit coupling. This limit is naturally incorporated in the 
theory proposed by Bolte and Keppeler ||] that based on the 7l — > 
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expansion in the Dirac (or Schrodinger) equation. Bolte and Keppeler 
have shown that the semiclassical trace formula without spin-orbit 
interaction acquires an additional modulation factor due to spin, but 
otherwise remains unchanged. We obtain the same modulation factor 
using the spin coherent state method. 



2 Classical dynamics and periodic 
orbits 

We begin with the classical phase space symbol of the Hamiltonian [|J 

H(p, q, z) = H Q (p, q) + KhSa(z) ■ C(p, q) = H + hH so . (1) 

It includes the spin-orbit interaction term fiH so which is linear in spin, 
but otherwise is an arbitrary function of (possibly multidimensional) 
momenta and coordinates p and q. The spin direction is described by 
a unit vector (t(z) = f (z\a\z), where <x are the Pauli matrices and 
the complex variable z = u — iv labels the spin coherent states of total 
spin S H. At the end of our calculations we will set S = \. The 
Planck constant appears explicitly in this classical Hamiltonian and is 
treated as the perturbation parameter in the weak-coupling limit. The 
spin-orbit coupling strength k is kept finite. Thus the condition fi — » 
provides both the semiclassical (high energy) and the weak-coupling 
limits. 

The Hamiltonian (Q) determines the classical equations of motion 
for the orbital and spin degrees of freedom H] 

OH dH dC 

P = — 7T- = — o KhSa-—, (2) 

oq oq oq 

dH dH 8C 

1 = = "5 hKhScr-— , (3) 

op op op 

a = -kct x C. (4) 

Since 

°"0) = 7-7^(>' 2u 'M 2 - 1 ) T ( 5 ) 
1 + \z\ A v ' 

in the "south" gauge,[] Eq. (^) is equivalent to two Hamilton-like equa- 



1 By the south gauge we mean the choice of parameterization of the spin coherent states 
by z such that <J z (\z\ — * oo) = 1. 
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tions 



(l + \z\ 2 ) 2 OH k , i2n2 da n 

In the leading order in fi we keep only the unperturbed terms in Eqs. 
@j and (H). It follows then that the orbital motion, in this approxi- 
mation, is unaffected by spin. The spin motion is determined by the 
unperturbed orbital motion via Eq. which does not contain h. It 
describes the spin precession in the time-dependent effective magnetic 
field C (po(t),qo(t)), where (po(t), qo(t)) is an orbit of the unperturbed 
Hamiltonian Hq. 

In order to apply a trace formula for the density of states, we need 
to know the periodic orbits of the system, both in orbital and spin 
phase space coordinates. The orbital part of a periodic trajectory is 
necessarily a periodic orbit of Hq. Assume that such an orbit with 
period To is given. Then Eq. (||) generates a map on the Bloch sphere 
<t(0) i — ► (t(Tq) between the initial and final points of a spin trajectory 
<r(t). The fixed points of this map correspond to periodic orbits with 
the period To. Since Eq. (||) is linear in <r, for any two trajectories 
(Ti(t) and (T2(t), their difference also satisfies this equation. But this 
means that \<Ti{t) — <r^{t)\ = const, i.e., the angles between the vectors 
do not change during the motion. Hence the map is a rotation by an 
angle a about some axis through the center of the Bloch sphere. The 
points of intersection of this axis and the sphere are the fixed points 
of the map (Fig. [I]). Thus for a given periodic orbit of Hq, there are 
two periodic orbits of H with opposite spin orientations (unless a is a 
multiple of 2tt, by accident). The angle a was mentioned in Ref. ||]. 



3 Modulation factor 
3.1 Correction to the action 

In order to derive a modulation factor in the trace formula, we need to 
determine the correction to the action and the stability determinant 
due to the spin-orbit interaction. The action along a periodic orbit 
is i 

/ n*- 1 u dv — vdu 

S = J pdq + 2SH <j> 1 + , , 2 = S pq + TuSspin. (8) 
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Figure 1: Axis of rotation and fixed points of the map <x(0) i — > <t(Tq 



While the spin part contains h explicitly, we need to extract the lead- 
ing order correction to the orbital action. This is the only place where 
we implicitly take into account the influence of spin on the orbital 
motion. It is convenient for the following calculation to parameterize 
both the perturbed and unperturbed orbits by a variable s S [0, 1], 



i.e., 



1 d ^ 

p — as. 

as 



(9) 



The time parameterization would be problematic since the periods of 
the perturbed and the unperturbed orbits differ by order of % (see 
Appendix [A]) . The correction to the orbital part due to the perturba- 
tion is 



SS 



' d 1° , d (X \ 
Op— +Po — {dq) 



o 



ds 



ds 



ds 



dq dp . 
op— °Q~T~ I " s +Po°Q 

Hi J CtO 



(10) 



The boundary term vanishes for the periodic orbit, and the integration 
can be done over the period of the unperturbed orbit now: 



5S,, 



To 



(5pq - 5qp ) dt 



To 



r 9H 
5p— — +. 
op 



dH c 
dq 



dt 



To 



5H dt. 



(11) 
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Since the perturbed and unperturbed orbits have the same energy, the 
variation of the Hamiltonian is 5Hq = —KH so . Taking into account 
Eq. (151), we can express the change in the orbital action as 



5S„ n = —hnS 



T 





( 2n \ 


C-adt = -HkS J ° C ■ 


2v 







dt 



(12) 

We now turn to the spin action. Parameterizing the trajectory 
with time and then using the equations of motion (|6|), (fa) and Eq. (|5|), 
we find 

TlkS [ t ° ( da da\ , , l2 , , 

T ° „ (*{}- \ Z Q \ dt 



UkS / C | v(l-\z\ 2 ) \ ——. (13) 



2\z\ 2 



1+ z 



Summing up the orbital and spin contributions Eqs. (12) and fll3|), we 
obtain the entire change in action as 

rT 



where 



5S = 5S pq + hS spill = hS [ ° F{t)dt, (14) 

Jo 



-u 

F(t) = kC • | -v | . (15) 



3.2 Stability determinant 

The stability determinant is derived from the second variation of the 
Hamiltonian about the periodic orbit |l|]. In the leading order 
in h, the orbital and spin degrees of freedom in are separated. 
This means that the spin phase space provides an additional block to 
the unperturbed monodromy matrix of the orbital phase space, which 
results in a separate stability determinant due to spin. The linearized 
momentum and coordinate for spin 



£ \ ( 5u 



v ) 1 + |z| 2 \ Sv 
satisfy the equations of motion 



(16) 



/ 



mi ; . (17) 
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Solving these equations we find the spin block of the monodromy ma- 
trix to be (Appendix ||) 



M 



where the stability angle is 



cos tp — sm tp 
sin (p cos <p 



T 



F{t)dt. 



(18) 



(19) 



The proportionality between <p and 5S [Eq. (|14])] will be exploited in 
a moment but, first, we find the stability determinant 



|det(M-/)| 1/2 = 2 
where / is the 2x2 unit matrix. 



. <P 
sm — 

2 



(20) 



3.3 Trace formula 

As was explained at the end of Sec. ^, for each unperturbed periodic 
orbit there are two new periodic orbits with opposite spin orientations 
<r(i). It is easy to deduce, then, that for these two orbits both 5S and 
(p have the same magnitude but opposite signs. Now we are ready to 
write the trace formula for the oscillatory part of the density of states 



An 



po ± 



■ COS 



1 7T 

-(So±SS)- -(mo + M±) 



(21) 



where the first sum is over the unperturbed periodic orbits and the 
second sum takes care of the contribution of the two spin orientations; 
Ao is the prefactor for the unperturbed orbit, which depends on the 
stability determinant and the primitive period; So and no are the 
unperturbed action and the Maslov index, respectively; \x± are the 
additional Maslov indices due to spin. The nature of spin requires 
an additional Kochetov-Solari phase correction Q that results in the 
shift S i — > S + \ of the total spin parameter in 5S (Appendix |C]). 
Setting 5 = ^, we end up with 

SiS — ► SS = hp. (22) 
With this relation and the additional Maslov index (Appendix |^) 



H± = 1 + 2 



[±-1 
L 2ttJ 



(23) 
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{[x] is the largest integer < x) the sum over the spin orientations in 
Eq. (|2l|) becomes 



— (,,,s 



± 2sin 2 



So 

n 



2 cos 



Aq cos 



7T 

So vr 




(24) 



This is our main result: each term in the periodic orbit sum is the con- 
tribution of an unperturbed orbit .4.0 cos — §/^q] times the modu- 
lation factor 



M = 2 cos 



(25) 



Note that no assumption was made on whether the unperturbed pe- 
riodic orbits are isolated or not. 



4 Comparison with another method 

Bolte and Keppeler ||] derived the modulation factor in the weak- 
coupling limit by a different method. Their results^ are expressed in 
terms of a spin trajectory with the initial condition 

er(O) = (0,0, -1) T (26) 

that obeys Eq. (Q). This trajectory, in general, is not periodic. As in 
our approach, the influence of spin on the orbital motion is neglected. 
The spin motion can be described by the polar angles (9{t), <fi(t)) with 
9(0) = ir. The modulation factor is then 

Mbk = 2 cos (^pj cos x, (27) 

where A9 = tt-9(T ) andP] 

X = -^J C-<rdt + -J [l + cos0(t)]0(t)cft. (28) 



In order to show that our modulation factor Eq. (|25| ) is equal to 
Mbk i let us express ip in terms of the polar angles. From Eq. (|5|) 
follows the coordinate transformation 

u = cot - cos q>, 
2 r ' 

Q 

v = cot - sin 4>. (29) 
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We reformulate the results of Ref. Q for the south gauge. 
3 Ref. (2| defines the phase 77 = — x- 
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Since tp oc 5S, we can represent it as a sum of two terms [cf. Eqs. 
©-©] 



/ C-adt+- [l + cosO(t)](j>(t)dt. (30) 
Jo * Jo 



2 ~ ~2 



There is a striking similarity between the expressions for x an d \ ■ The 
only difference is that in Eq. (j28|) the integration is, in general, along a 
non-periodic orbit with the initial condition Eq. fl26|) , while in Eq. fl30| ) 
the integration is along the periodic orbit. Since the modulation factor 
should not depend on the choice of the z direction, we can choose the 
z axis to coincide with the spin vector <r(0) for the periodic orbit at 
t = 0, i.e., the z axis will be the rotation axis in Fig. [|. Then one of the 
periodic orbits will satisfy the initial condition Eq. fl2q), and thus both 
X and ^ can be calculated along this orbit and are equal. Moreover, 
A9 = in this case. Therefore the modulation factors derived within 
the two approaches coincide, 

M B k=M. (31) 

It was mentioned in Ref . S that Ai bk = 2 cos % , where a is the 
rotation angle defined in Sec. |2[ Then, of course, we conclude that 

cos ^ = cos ^. (32) 



To see that this is indeed the case, we can go back to Sec. 3.2 where we 
calculated the stability determinant. It follows from that calculation 
that the neighborhood of the periodic orbit is rotated by an angle ip 
during the period (Appendix |B|). Therefore the entire Bloch sphere is 
rotated by this angle. Clearly, the angle of rotation must be defined 
mod 4-7T, i.e., it depends on the parity of the number of full revolutions 
of the Bloch sphere around the periodic orbit during the period. It 
would be desirable to prove Eq. ( |3~2"|) without referring to the small 
neighborhood of the periodic orbit. 

The same property can be also shown if one treats the spin quan- 
tum mechanically. The spin propagator for the choice of the z axis 
along the rotation axis (so that % = 2 ) 1S @ 



rf(To) = ( 6 o* Jf ' (33> 



Applying this operator to a spinor 

(V^,^-) at i = 0, we get the 
spinor I ip + e~ l 2 , tp^e 1 ^ I at t = Tq, which corresponds to the initial 



spin vector rotated by the angle ip about the z axis, i.e., 99 = a. 
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5 Summary and conclusions 



We have studied the case of weak spin-orbit coupling in the semiclas- 
sical approximation using the spin coherent state method. The limit 
is achieved formally by setting Ti — > 0. The trajectories in the orbital 
subspace of the extended phase space then remain unchanged by the 
spin-orbit interaction. For each periodic orbit in the orbital subspace 
there are two periodic orbits in the full phase space with opposite 
spin orientations. The semiclassical trace formula can be expressed 
as a sum over unperturbed periodic orbits with a modulation factor. 
This agrees with the results of Bolte and Keppeler. The form of the 
modulation factor does not depend on whether the unperturbed sys- 
tem has isolated orbits or whether it contains families of degenerate 
orbits due to continuous symmetries. 

We remark that in the semiclassical treatment of pure spin sys- 
tems, a renormalization procedure is necessary in order to correct the 
stationary phase approximation in the path integral for finite spin S. 
Such a renormalization is equivalent to the Kochetov-Solari phase cor- 
rection that we employed here without justification for a system with 
spin-orbit interaction. Although this correction worked well in our 
case, it may be necessary to develop a general renormalization scheme 
when the interaction is not weak. 
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A Time parameterization 



For pedagogical reasons we do the calculation in Eq. (10) with the 
time parameterization. In this case S pq = J pqdt, where T is the 
exact period. Then the correction is 



p<i 



f T ° r 
Jo i' 



T 



Sp% + Po(Sq) dt + p (T )q (T )5T 

To 



(5pq - 5qp )dt + p c 



+ p (T )q (T )5T. (34) 
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Transforming the boundary term 



To 

= p (T ) [q(T ) - q (T ) - q(0) + q (0)} = p (T ) [q(T ) - q(0)} 



= po(T ) [q(T ) - q(T)) ~ -p (T )q (T )5T, (35) 



we see that it cancels the period correction term. 



B Monodromy matrix 



We derive the monodromy matrix Eq. (|18[) . In order to solve the 
equations of motion ( |T7| ) we define r\ = £ + iv. Then rj = ir)F(t), 
which solves to 



rj(t) = 77(0) exp 

It follows then that 



t 



% \ F(t')di! 







(36) 



£(Tq) = ^(0) cos if — i/(0) sin 99, 

i/(T ) = ^(0) sin if + i/(0) cos <p, (37) 



resulting in Eq. (JX8J) . 

Note that according to Eq. 



l-cr 2 



1/ = V^S (5a v + . (38) 

If we choose the z axis in such a way that the periodic orbit starts 
and ends in the south pole, i.e., <x(0) = <t{Tq) = (0, 0, — 1) T , then at 
t = and t = To we have 

£ = VhS5a x , 

v = VhS5a y . (39) 



Comparing with Eqs. (37) we conclude that the neighborhood of the 



periodic orbit is rotated by the angle tp after the period. 



C Kochetov-Solari phase shift 

The Kochetov-Solari phase shift M is given by 



1 



T 



ip KS = - / A(t)dt, (40) 
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where 



A(t) 



1 

2h 







[l + \z\ 2 ?dH 



dz 



2S 



dz 



+ c.c. 



(41) 



The spin-dependent part of the Hamiltonian is [cf. Eq. 



HH so (z,z) 



c 



i\z 



z 

■z) 



(42) 



We find 



0(1 + 



dz 



\z\ 2 ) 2 dH« 



2S 



dz 



kC- 



—z 
iz 
1 



(43) 



therefore the phase shift becomes 



(44) 



Comparing to Eq. (14) we see that it effectively shifts the spin S by 
|. One should keep in mind that this phase correction was originally 
derived for a pure spin system. It has not been proven to have the 
same form for a system with spin-orbit interaction. In the special case 
of the weak-coupling limit we have a reason to believe that the result 
Eq. ( [44|) is correct, since we were able to reproduce the modulation 
factor found with another method (see Sec. ||). 



D Maslov indices 

The additional Maslov indices \x± are determined by the linearized 
spin motion about the periodic orbit. The second variation of the 
Hamiltonian reads [cf. Eq. fll^)] 



H^(^) = ^ 1 (e + v 2 )- (45) 



Following Sugita || we define its normal form 

Hnorm = ^{e + V 2 ) (46) 

that has a constant frequency and generates the same phase change 
if as after the period To. Then the spin block of the monodromy 
matrix can be classified as elliptic and its Maslov index is given by 
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Eq. (^). 93 is the stability angle of one of the two orbits with opposite 
spin orientations. Therefore, without loss of generality, we can assume 
that ip > 0. Then, explicitly, 

f ±1, if v3G(0,2tt) mod47r 
M± \ ±3, if ip£{2tt,4tt) mod 4^ ' 1 j 

On the other hand, 

sign (sin £ W if ^ G ( °' 27r) m ° d ^ (48) 

S g V S 2 J if (27r,47r) mod 4^ ' [ ^ } 

Clearly, one can take fi± = ±1 and at the same time remove the 
absolute value sign from sin ^ , as was done in Eq. . 
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